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REVISION 3 CHANGES 


The following is a summary of the major technical changes included in this 
revision: 

1. A plane change maneuver has been incorporated in 
the program's logic. 

2. In the search for an apsidal crossing, program checks 
have been included to update through intervals of v 

in the case of near circular orbits or close proximity 
to an apsidal point. 

3. Following an astronaut overwrite of a Lambert maneuver, 
the state vector is updated to establish a new target 
vector for use in the powered flight routine. 

4. Following the computation of each maneuver, the 
position vector of the primary vehicle is offset — 
to help compensate for the effects of the finite ma- 
neuver - prior to the update of the state vector to the 
next maneuver point. 
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FOREWORD 


This document is one of a series of candidates for inclusion in a future re- 
vision of JSC-04217, "Space Shuttle Guidance, Navigation and Control Design Equa- 
tions. " The enclosed has been prepared under NAS9- 10268, Task No. 15- A, 

"GN & C Flight Equation Specification Support", and applies to functions 1, 2, 3, 
and 6 of the Rendezvous Targeting Module (OG3) as defined in JSC -03690, Rev. D, 
"Space Shuttle Orbiter Guidance, Navigation and Control Software Functional Require- 
ments", dated January 1973. 
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NOMENCLATURE 


Semi-major axis of a conic 
Alarm code i 

Estimated magnitude of thrust acceleration 

Failure in fuel optimization loop 

Failure in height loop 

Failure in phasing loop 

Failure in obtaining Lambert solution in 
General Maneuver Routine 

Failure to find perigee/ apogee in Search Routine 

Incompatible altitudes and elevation angle 

Failure to find time corresponding to elevation 
angle in Search Routine 

Failure to find desired position vector in 
Desired Position Routine 

Failure to update through 0 in Update Routine 

Iteration counter 

Height iteration counter 

Phase iteration counter 

Intermediate variables 

Delta altitude 


Delta projected position 



At 

At 


max 


Av 


A -LOS 


A — LV 


A v. 


Av 


A v 


T 


Ax 


a e 


Delta time 

Maximum time step allowed in Search Routine iteration 
Maneuver velocity 

Maneuver in line-of-sight coordinates 

Maneuver in local vertical coordinates 

Av used during height maneuver 

Av used during phasing maneuver 

Delta velocity used in fuel minimization loop 

Delta independent variable 

Delta central angle 


e 


e h 

e p 

e L 

§0 

i 

-N 


i 

1 max 


m 

M 

n 


Error 

Eccentricity 
Height error 
Phasing error 

Elevation angle (defined in Figure 2 ) 

Reference value for the acceleration of gravity 
Unit vector 

Unit normal to the plane used in powered flight 
guidance 

Number of the maneuver 

Number of the last maneuver in the 
maneuver sequence 

Estimated vehicle mass 

Rotational matrix 

Vector normal to the orbital plane 

Number of revolutions 


vii 



p 

Partial used in Newton Raphson 


iteration 

r 

Distance ratio 

r 

Position vector 

— D 

Desired position vector 

lie 

Target vector used in powered flight guidance 

s 

Switch used in Desired Position Routine 

s 

astro 

Astronaut overwrite switch 

S coplan 

Coplanar switch 

s direct 

Maneuver direction switch 

S eng 

Engine select switch 

s exit 

Program exit switch 

5 fail 

Failure switch 

s man 

Maneuver switch 

S opt 

Maneuver optimizing switch 

s outp 

Out -of -plane switch 

S pert 

Perturbation switch 

3 phase 

Phase match switch 

S proj 

Projection switch 

s rdes 

Desired position switch 

s soln 

Solution switch 

s search 

Search switch 

S tar 

Target offset switch 

3 term 

Terminal constraint switch 

s update 

Update switch 
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Time 


Final time 
Velocity vector 

Vertical component of velocity 
Circular velocity 

Independent variable in Iteration Routine 
Out-of-plane parameters (see Figure 6a) 

Radial component of velocity divided by v^ 

Horizontal component of velocity divided 
by v c 

Tolerance on time in fuel optimizing loop 

Tolerance on height in height loop 

Tolerance on central angle in phasing loop 

Tolerance on central angle in Search Routine 

Tolerance on elevation angle in Search Routine 

Altitude increment in Search Routine 

Tolerance on central angle in Desired Position Routine 

Tolerance on central angle in Update Routine 

Flight path angle 

Gravitational constant 

Central angle 


Perigee angle 



Subscripts 


F 

i 

I 

LOS 

LV 

N 

0 

P 

s 

T 

TA 


Final 

Number of the maneuver 
Initial 

Line -of -sight 
Local vertical 
New 
Old 

Primary 

Stored 

Target 

Target for primary vehicle 


x 
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INTRODUCTION 


The rendezvous of the Shuttle (primary vehicle) with a target vehicle (e. g. 
the Space Station) is accomplished by maneuvering the Shuttle into a trajectory 
that intercepts the target vehicle orbit at a time that results in the rendezvous of 
the two vehicles. The function of rendezvous targeting is to determine the targeting 
parameters for the powered flight guidance for each of the maneuvers made by the 
Shuttle during the rendezvous sequence. 

In order to construct the multimaneuver rendezvous trajectory, sufficient 
constaints must be imposed to determine the desired trajectory. Constraints 
associated with the Shuttle mission will involve such considerations as fuel, light- 
ing, navigation, communication, time, and altitude. The function of premission 
analysis is to convert these — which are generally qualitative constraints-into a set 
of secondary quantitative constraints that can be used by the onboard targeting 
program. By judicious selection of the secondary constraints, it should be possible 
to determine off-nominal trajectories that come close to satisfying the primary 
constraints. 

The proposed onboard rendezvous targeting program consists primarily ol a 
main program and a generalized multiple-option maneuver subroutine. The driving 
program automatically and sequentially calls the maneuver subroutine to construct 
the rendezvous configuration from a series of maneuver segments. * The main pro- 
gram is capable of handling rendezvous sequences involving any given number of 
maneuvers. Enough different types of maneuver constraints are incorporated into 
the subroutine to provide the flexibility required to select the best set of secondary 
constraints during premission planning. In addition, the astronaut has a large, 
well defined list of maneuver options if he chooses to modify the selected nominal 
rendezvous scheme. 

As the new approach represents, in essence, just one targeting program, 
there is considerable savings in computer-storage requirements compared to former 
approaches in which each maneuver used in the rendezvous scheme had a separate 
targeting program. The programming and verification processes of this unified 
approach will also result in implementation efficiencies. 

1 . 1 Number of Independent Constraints Involved in a 

Rendezvous Sequence 

During the Gemini and Apollo flights and in the design of the Skylab rendezvous 
scheme various numbers of maneuvers were utilized in the rendezvous sequence. 

The range went from two (Apollo 14 and 15) to six (Skylab). 
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The number of independent constraints (i.e., the number of explicitly satisfied con- 
straints) in each rendezvous sequence must equal the number of degrees of freedom 
implicitly contained in the sequence. To establish this number, a rendezvous con- 
figuration can be constructed by imposing arbitrary constraints until the configura- 
tion is uniquely defined. For example, a four maneuver coplanar sequence is 
shown in Figure 1, followed by a coast to a terminal point. Using the constraints 
V- (velocity magnitude), r^ and 0^, it is easy to establish that the total number 
involved is 12, assuming the time of the first maneuver has been established. Re- 
moving one maneuver will reduce the number of degrees of freedom by three. 

Hence, the number of independent constraints necessary to uniquely determine the 
maneuver sequences are 


Number of maneuvers 
in sequence 

1 

2 

3 

4 
etc 


Number of independent 
constraints required 

3 

6 

9 

12 


If the above rendezvous are not coplanar, one additional constraint has to be added 
to each sequence to allow for the out-of-plane component. 

In some cases the number of primary constraints may be insufficient to 
uniquely determine a rendezvous trajectory for the desired number of maneuvers. 
One way of overcoming this deficiency in constraints is by introducing sufficient 
variables to complete the determination of the rendezvous trajectory and then 
determining values for these variables by minimizing the fuel used. 

In order to take advantage of updated state vectors due to navigation or ground 
updates, the rendezvous targeting program is called prior to each maneuver to 
compute the upcoming maneuver. In general, each maneuver computation will in- 
volve a multimaneuver sequence as the nature of the targeting constraints do not 
allow the maneuvers to be independently computed." The relationship between the 
rendezvous sequence involving n maneuvers and the maneuver sequences is shown 
below. Maneuver Segments 


Rendezvous Sequence 


Maneuver Sequences 


1 2 3 4 n- 1 



(maneuver points) 
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Each maneuver sequence is composed of a number of maneuver segments and is 
basically independent from the other maneuver sequences. These sequences must 
have the same number of independent constraints as tabulated above. 

1. 2 The Construction of a Maneuver Segment 

Each n-maneuver rendezvous sequence can be divided into n-1 maneuver seg- 
ments. Each segment involves, basically, the addition of a maneuver to the primary 
vehicle's velocity vector and an update of both vehicle's state vectors to the next 
maneuver point. 

A maneuver segment is herein generated in one of three ways: 

Forward generation A maneuver A v is computed and added 

to the velocity vector in a specified di- 
rection. The state vector of the primary 
vehicle is then updated through a speci- 
fied amount to arrive at the next maneu- 
ver position. 

Target generation The target vehicle is updated through a 

specified amount and then offset to establish 
a target vector. An option is available at 
this point to compute a coelliptic velocity 
vector and update through A t to establish 
a new target vector as shown below. 



The maneuver is then computed by uniquely 
specifying the nature of the traverse between 
the primary vehicle's position and the target 
vector. 
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Integrated generation 


In this case, the maneuver segment 
is computed as an integral part of a 
maneuver sequence involving more 
than one maneuver segment. The 
nature of the constraints are such that 
the maneuver sequence cannot be sub- 
divided into uniquely defined maneuver 
segments. The maneuver segment will 
usually have one degree of freedom, 
which will generally be assumed to be 
the magnitude of the maneuver. 

Each of the above methods is defined by specifying trajectory constraints 

by setting certain switches and parameter values. Specifying a trajectory constraint 

is equivalent to specifying one or more independent constraints. On the other hand, 

specifying an independent constraint can also be equivalent to specifying one or more 

trajectory variables. (See Hef. 9) A trajectory constraint common to all three of 

the above methods is the state vector update switch s . . The options associated 

update 

with this switch are: 

Update from time t to time t,-, 

r 

Update through time interval At 

Update through n revolutions 

Update through 0 radians 

Update to be colinear with a 
specified position vector 


s = s 

update 


In the next three sections, the trajectory constraints associated with each of the 
above methods will be listed. 
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1.2.1 Maneuver Options in Forward Generation of 
Maneuver Segment 


The forward generation of a maneuver segment is accomplished in one of two 
ways. Either the maneuver magnitude is uniquely determined in terms of the state 
vector at the maneuver time or the maneuver is determined by an iterative search 
to satisfy a terminal constraint. 

The maneuver magnitude Av is either calculated or assumed depending on 

the maneuver switch s , and it is applied in a direction controlled by the di- 

man 

rection switch s .. The options associated with the maneuver switch are: 


1 

2 


s 

man 





3 

4 


A v is assumed specified 

Av is computed based on a post maneuver 
velocity vector being "coelliptic" with the 
state vector of the target vehicle 

A v is computed from the conic circular 
velocity constraint 

A v is computed based on a Hohmann type 
transfer resulting in a Ah change in 
altitude 


The options associated with the maneuver direction switch are : 


-1 


s direct 




-2 


Apply Av is horizontal direction in plane 
of primary vehicle 

x^pply Av in horizontal direction parallel 
to orbital plane of the target vehicle 

Apply Av along velocity vector in plane 
of primary vehicle 

Apply Av along velocity vector parallel 
to orbital plane of the target vehicle 


The selection of the update switch s update determines the update of the 
primary vehicle's trajectory following the maneuver to the position of the next 
maneuver. A terminal constraint can be imposed at this point by setting the ter- 
minal switch s : 
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term 



Terminal constraint is a height constraint 
Terminal constraint is a phasing constraint 


Following the computation of the height/ phasing error, the maneuver magnitude is 
varied in an iterative search to satisfy the height/ phasing constraint. 


1.2.2 Maneuver Options in Target Generation 
of Maneuver Segment ~ 

The target generation of a maneuver segment starts with the selection of 
the update switch for the target vehicle. If this switch equals four, 0 will be 
augmented by the central angle between the primary and target vehicles before 
being used. The position of the target vehicle is then offset through either (e^, Ah) 
or (AS, Ah), depending on whether s^, is negative or positive, to obtain a target 
vector. The "TPI offsets" (e^, Ah) are discussed in Section 5 (see Figure 2 for 


definition of e^). 


If 


s ^ | equals two, a coelliptic velocity vector is computed 
tcLI* 


based on the target vector, and a new target vector is defined by updating the co- 
elliptic state vector through At. The options associated with s tar are: 


-2 Offset target (e L> A h ). Compute coelliptic 
velocity and update through (negative) At . 

-1 Offset target (e^, Ah) 

s tar = 

0 No target offset 

1 Offset target { A0 , Ah) 

2 Offset target (A0, Ah). Compute coelliptic 

• [_ velocity and update through (negative) At 

The nature of the traverse between the primary 

vehicle's initial state vector and the target vector is 

controlled by the maneuver switch s : 

J man 

5 Lambert - the trajectory is time constrained 

6 Horizontal - the maneuver is constrained to 
be in the horizontal direction 

s = « 

7 Tangential - the maneuver is constrained to 
be in the direction of the velocity vector 

8 Apogee/Perigee - the trajectory has an apogee/ 
perigee occurring at the target point. 

There is a minimum Av option associated with the above maneuvers which 
is controlled with the optimization switch • 
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Position Vector of Primary Vehicle 


[ = Unit horizontal in forward direction for primary vehicle 

LOS = Line of Sight 

1. If the LOS projection on Ms positive*. 

a. When the LOS is above the horizontal plane. 0 < eL < 7r/ 2 

b. When the LOS is below the horizontal plane, 37r/ 2 < e L < 2 tt 

2. If the LOS projection on [ is negative.- 

a. When the LOS is above the horizontal plane, irl 2<ej_ < zr 

b. When the LOS is below the horizontal plane, n<ei < 37 r/ 2 


Figure 2. Definition of the Elevation Angle e^ 
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-2 Minimize the sum of the magnitude of 
the first and the next maneuver (based 
on a coelliptic parting velocity) by 
varying At , the time of update of the 
target vehicle. 

-1 Minimize the magnitude of the first ma- 
neuver by varying £t, the time of update 
of the target vehicle. 

s = 

^ 1 Minimize the magnitude of the first ma- 

neuver by varying At, the time between 
the next maneuver and the initial offset 
position. (See sketch on page 1-4 ) 

2 Minimize the sum of the magnitudes of 
the first and the next maneuver (based 
on a coelliptic parting velocity) by 
varying At, the time between the next 
maneuver and the initial offset position 
(see sketch on page 1-4). 

This minimization is accomplished by driving the slope (Av / independent variable) 
to zero using a Newton Raphson iteration scheme. 

i-2.3 Maneuver Options in Integrated Generation 
of Maneuver Segment 

The integrated generation of a maneuver segment involves an iterative solu- 
tion to determine a maneuver sequence which cannot be sequentially solved for its 
maneuver segment components. The maneuver is computed by guessing its magni- 
tude, assigning a direction and plane through selection of the direction switch 
S direct' u P datin S the primary vehicle’ s state vector after selecting switch s 

update 

and then calling additional maneuver segments until reaching the point at which P 
the terminal constraint is to be attained. The maneuver is then iteratively de- 
termined by satisfying the terminal constraint. The number of additional maneu- 
ver segments and the nature of the terminal constraint are controlled by the ter- 
minal constraint switch s 

term 
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The terminal constraint is a phasing 

constraint and it occurs at the I s I 

| term[ 

maneuver point from the start of the 
maneuver segment . 

The terminal constraint is a height 
constraint and it occurs at the s 

term 

maneuver point from the start of the 
maneuver segment. 

Both a height and phasing constraint 
occur at the same maneuver point. 

The first digit m of s^ repre- 
sents a phasing constraint that occurs 
at the n ^ maneuver point from the 
start of the phasing maneuver segment. 
The last digit n 2 of s term represents 
a height constraint that occurs at the 
n 2 maneuver point from the start of 
the height maneuver segment. 


1.2.4 Summary of the Maneuver Segment Con straints 

The maneuver and trajectory constraints that can be imposed on a maneuver 

segment can be divided into the following catagories (see Figure 3). 

Primary vehicle update constraints 
Target vehicle update constraints 
Initial velocity constraints 
Offset constraints 
Terminal constraints 
Traverse constraints 


Table 1 contains a detailed listing of the constraints. The three independent con- 
straints (four in the case of noncoplanar traverses) which govern a maneuver seg 
ment cannot be chosen arbitrarily from this list for two reasons: 

(1) There is not a one-to-one correspondence. between the 
trajectory constraints and the independent constraints. 
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(2) Selecting some constraints negates the need for some 
others (e.g. selecting a Lambert constraint negates the 
need for a maneuver direction constraint). 

In the case of a straight forward rendezvous profile, a basic understanding of the 
nature of the constraints should allow the constructor of the rendezvous sequence 
to choose a set of trajectory constraints which determine the required number of 
independent constraints. For a complex rendezvous profile, such as Skylab a 
more formal approach such as presented in Reference 11 should be used. One of 
the justifications for presenting the three methods of generating a maneuver seg- 
ment was to aid the constructor of the rendezvous sequence in choosing compatible 
sets of constraints. 
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TABLE 1 


DETAILED LISTING OF CONSTRAINTS 
(Sheet 1 of 2) 

Primary and Target Vehicle Update Constraints 
Delta time 
Initial and final time 
Central angle 
Number of revolutions 
Terminal position vector 

Initial Velocity Constraints 
Plane 

Parallel to target orbit 
Parallel to primary orbit 

Direction 

Horizontal 

Along velocity vector 
Magnitude 

Circular 
Coelliptic 
Altitude change 
Specified 

Offset Constraints 

Angle 
Altitude 
Elevation angle 

Terminal Constraints 
Height 
Phase 
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TABLE 1 


DETAILED LISTING OF CONSTRAINTS 
(Sheet 2 of 2 ) 


Traverse Constraints 
Minimum Fuel 

One maneuver optimization 
Two maneuver optimization 
Apogee/Perigee designation 
Horizontal maneuver 
Tangential maneuver 
Lambert (time) 


1-14 



2 . 


FUNCTIONAL FLOW DIAGRAMS 


The rendezvous targeting program consists of two major parts — a generalized 
maneuver subroutine which basically computes a maneuver and updates the state 
vectors of both vehicles to the time of the next maneuver and a main program which 
sequentially calls the subroutine to assemble a rendezvous sequence. These pro- 
grams call a number of subroutines which are briefly described below and in de- 
tail in Section 5. 


Search 


Phase 

Match 


Desired 

Position 


To update the state vectors to either a 
specified apsidal crossing, a time, or 
an elevation angle. 

To phase match the target vehicle's 
state vector to the primary vehicle' s 
position vector. 

To compute an offset target vector or a desired 
position to be used in a phasing constraint. 


Update 


To update a state vector through a speci 
fied interval. 


Coelliptic - To compute a coelliptic velocity vector. 
Maneuver 


Iteration - To determine a new estimate of the in- 
dependent variable in a Newton Raphson 
iteration scheme. 


The functional flow diagram for the main program is shown in Figure 4. The 
main function of this program is to sequentially call the General Maneuver Routine 
to compute each maneuver segment for maneuvers numbered from i to i- max * 
There are three major options that can be exercised prior to the calculation of the 
first maneuver segment: 

(1) A search for the time of the first maneuver. 

This time can be specified by: 

(a) An elevation angle, which is to be 
attained at the maneuver time. 
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Figure 4 . Main Program - Function Flow Diagram 
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(b) Whether the next maneuver should occur at 
the next apsidal crossing, the next perigee 
crossing or the nth apsidal crossing. 

(2) A phase matching of the state vector. 

(3) A rotation of the primary vehicle' s state vector 
into the plane of the target vehicle. 

There are three separate iterative loops built around the call to the general 
maneuver routine. One loop serves to minimize the fuel used during a maneuver 
segment with the options determined by the optimizing switch. 

The other two iterative loops involve maneuver segments which contain con- 
straints that do not allow the explicit calculation of the maneuver. These con- 
straints are height and phasing constraints imposed at the end of a maneuver seg- 
ment and controlled with the terminal switch. The iterative loop will involve 
several maneuver segments if sufficient constraints are not imposed to solve each 
segment uniquely. 

The functional flow diagram for the general maneuver routine is shown in 
Figure 5. This routine generates the departure velocity at the initial point in one 
of two ways: 

(1) As an explicit function of the initial state vectors. 

(2) By defining a target vector and then computing an 
intercept trajectory based on a specified constraint 
(as indicated by the setting of S man >. The target 
vector is determined by offsetting the updated posi- 
tion vector of the target vehicle. Depending on the 
setting of the switch s tar , a coelliptic velocity vec- 
tor is computed at the offset point and the coelliptic 
state vector is updated through At to obtain a target 
vector. 

Following an update of both vehicle's state vectors to the time of the next maneu- 
ver, the Av used or the terminal height/ phase errors are calculated as required. 


2-3 



ENTER 



Figure 5. General Maneuver Routine - Functional Flow Diagram 
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3. 


INPUT AND OUTPUT VARIABLES 


The inputs to the orbiter rendezvous targeting program can be divided into 
five catagories. 

Pre-Maneuver Switches 

Upon selecting a maneuver from the rendezvous sequence, these switches 
(specified for each maneuver) serve in determining the state vectors at the maneu- 
ver point, the out-of- plane parameters and the calculation of a desired position 
vector. These inputs can also be used in determining the time of a specified apsidal 
crossing or the time at which a specified elevation angle is to be attained. 


s 


coplan 



s exit 



s outp 



S pert 



r 


0 

1 


’phase 



-2 


Coplanar switch 
Bypass 

Rotate primary state vector into plane 
of target vehicle's orbit 

Exit switch 

Bypass 

Exit from routine 
Out- of- plane switch 
Bypass 

Compute out-of-plane parameters 
Compute out-of-plane parameters and 
modify maneuver by -y p 

Perturbation switch 
Do conic state vector updates 
Include oblateness based on J 2 
Other perturbations as required 

Phase match switch 
Bypass 

Phase match state vectors (target leading 
primary) 

Phase match state vectors based on target 
leading primary by more than 360° 

Phase match state vectors (primary 
leading target) 

Phase match state vectors based on primary 
leading target by more than 360° 
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Desired position switch 



Bypass 

Compute desired position vector based on (e^, Ah) 
Compute desired position vector based on (A0, Ah) 


Search switch 


-4 

-3 

s search \ ^ 

-1 

0 

n 


Compute elevation angle 

Search for elevation angle 

Update to time L 

Search for next perigee crossing 

Bypass 

Search for the nth apsidal crossing 
(n>0) 


Maneuver Switches 

These switches (specified for each maneuver) set the constraints employed 
in determining the maneuver segments. 


s direct \ 


-2 

-1 

0 

1 


v. 


s 


man 



Direction switch 

Av in direction of primary's velocity 
vector, parallel to primary's orbital 
plane 

Av in horizontal direction, parallel 
to primary's orbital plane 
Bypass 

Av in horizontal direction, parallel to 
target's orbital plane 
Av in direction of primary's velocity 
vector, parallel to target's orbital plane 

Maneuver switch 

Av is specified 

Av is based on coelliptic velocity 
Av is based on circular velocity 
Av is based on altitude change 
Lambert maneuver to offset target vector 
Horizontal maneuver to offset target vector 
Tangential maneuver to offset target vector 
Perigee/ apogee insertion at offset target vector 
Plane change 
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Maneuver optimizing switch 
0 Bypass 

\ 1 Minimize Av^ 

2 Minimize Av. + Av. , , 
v 1 i + 1 


s 


soln 


Multi-revolution solution switch 

[-1 Solution with smallest initial flight path 
J angle (measured from local vertical) 

j 1 Solution with largest initial flight path angle 


' -2 Offset target (e T , Ah). Compute coelliptic 

i-J 

velocity and update through (negative) At 
. -1 Offset target (e^, Ah) 

0 No target offset 

1 Offset target ( A9 , Ah) 

2 Offset target ( A9 , Ah). Compute coelliptic 
velocity and update through (negative) At 


’term J 




Terminal constraint switch 

(<0) Compute phasing error and back up 
-(n + 1) maneuvers for start of phase loop 
Bypass 

(0 < n<10) Compute height error and back up 
n - 1 maneuvers for start of height loop 
(10 < n < 100) Phase and height loop terminate 
on same maneuver. For phase loop back up 
x - 1 (where x is first digit of n) maneuvers 
for start of phase loop. For height loop back 
up y - 1 (where y is last digit of n) maneuvers 
for start of height loop 


Update switch 


s 


update 



Update through tp - t 

Update through At 

Update through n r 

Update through 0 

Update to be colinear with r D 


Parameter Values 

The parameter values (specified for each maneuver) are values for the con- 
strained parameters. 
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Ah 

Delta altitude 

A hp. 

Delta altitude, final 

A9 

Delta central angle 

At 

Delta time 

A v 

Maneuver magnitude 

n r 

Number of revolutions 

*F 

Final time 

e L 

Elevation angle 


Post-Maneuver Switch 

This switch determines the options available following the calculation of the 
maneuver. 


s - \ 
astro ' 


Astronaut overwrite switch 
Bypass 

Overwrite maneuver in local vertical 
coordinates 

Overwrite maneuver in line -of- sight 
coordinates 


Maneuver Call Variables 


The maneuver call variables have to be specified for each call to the maneu- 
ver sequence. 

State vector of the primary vehicle 
State vector of the target vehicle 
Maneuver number 


— P J -P 


— T ’ -T 


t . 

1 


eng 

m 


Current time 


.th 


Time of the i maneuver (if s search J -2, then t. is not specified) 
Engine select switch 
Estimated vehicle mass 


Depending on the rendezvous sequence, there may also be some switches that have 
to be modified as a function of the maneuver number. 

Excluding the maneuver call variables, all the input variables can be set 
prior to the flight. 

The output parameters for the initial maneuver in the sequence are more 
complete than for the succeeding maneuvers. 
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Output Parameters for the Initial Maneuver 


Av. 

1 

A -LOSi 

A ^LVi 


-lc 


i 


N 


t 


T 


Maneuver magnitude 

Maneuver in line of sight coordinates 

Maneuver in local vertical coordinates 

Target vector used in Powered Flight Guidance 
Routine (See Ref. 5) 

Unit normal to plane used in same routine 
Time associated with 


Other parameters such as delta altitude, phasing angle, elevation angle and perigee 
altitude can be computed as required. 

Output Parameters for the Other Maneuvers in the Sequence 
t Time of the maneuver 

Av Maneuver magnitude 

Illustration of Inputs 

Table 2 contains a set of inputs for the Shuttle targeting program based on the 
five maneuver Skylab rendezvous configuration. The following switches and para- 
meters are not used as inputs to the Orbiter program: 

s astro 3 s exit ’ s opt’ s outp' s soln ' 3 


l 

The inputs in Table 2 are set prior to the mission so they will not have to be in- 
serted by the astronaut. The astronaut will have to modify the following quantities 
upon resetting the maneuver number as well as inserting the time of the next ma- 
neuver. 


i = 2: 

s, =0, 

term 2 

s , =32 

term^ 

i = 3: 

s =5, 

man^ 

At 3 = ' At NSR-TPI 

i = 4: 

s , =0 

term 4 
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TABLE 2 


INPUT VARIABLES FOR SKYLAB RENDEZVOUS CONFIGURATION 


Maneuver 


Input 

Variable 

1 

(NCI) 

2 

(NC2) 

' 3 

(NCC) 

4 

(NSR) 

5 

(TPI) 

s 

1 

1 

0 

1 

0 

coplan 






a 

1 

1 

1 

0 

0 

direct 






s 

1 

1 

1 

2 

5 

man 






s 

1 

1 

I 

0 

1 

pert 






s 

1 

• 1 

0 

0 

0 

phase 






Q 

-1 

-1 

0 

0 

0 

rdes 






g 

-2 

-2 

-2 

-2 

-3 

search 






s 



-2 


0 

b tar 






s 

0 

1 

0 

42 

0 

term • 






s 

3 

3 

2 

1 

2 

update p 






c 

update^ 



1 


4 

e L 

e L 

e L 

e L 


e L 

i 

4 

4 

3 

4 

5 

max 



* 



n 

r 

n rNC 1-NC2 

n rNC2-NCC 



V 

t TPI 

t TPI 

t TPI 

t TPI 


0 





9 TPF-TP1 

A t 



At NSR-NCC 



Ah 


Ah NCC 

Ah TPI 

Ah TPI 


Ah F 

A^tpi 

Ah Tpi 




Av 

Av NC1 

Av NC2 

Av ncc 
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4. 


DESCRIPTION OF EQUATIONS 


The only equations contained in this document which are not trivial are those 
involved in computing the traverse between two specified position vectors. The 
required equations can be derived from the equation of the conic expressed in the 
form 

r = r F / r i = 0j 2 / + e c cos ( 6 + 6 p ) J 

where 

22 2 2 
e c = l «i Bf + ( 0£ - i y ) 

= cos" 1 [( jS \ - 1 )/ e. c j (perigee angle) 
v c = W rj ) 1 / 2 

Oij and pj are the normalized (with respect to v c ) radial and horizontal components 
of velocity. 

The above equation can be expressed 

p s / r i r ~ sin 9 ~ c 2 ) / C 1 (1) 

where 

Cj = cos G - l/r 

c 2 = 1 - cos 6 

p = semilatus rectum 

For a maneuver that is constrained to be in a horizontal direction, Eq. (1) 
can be solved for Pj as a function of the specified . 

^1= [ o-j. sin 0 + (dj 2 sin 2 G - 4 c 2 ) 1/ ^ 2 ] / 2 

As there has to be both a positive and negative ^ solution to this equation (one 
trajectory in each rotational direction), the sign choice is resolved in favor of plus 

For a maneuver that is applied along the velocity vector, the flight path angle 
y j. is to be held fixed. Using Eq. (1) 
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Therefore, 


jSj = [ c 2 /(sin 8 tan y l - Cj)] 1 ^ 2 

By interchanging the I and F subscripts, Eq. (1) can be expressed 
p g = rptepgpSin 8 - c 2 )/ (cos 8 - r ) 

Combining with Eq. (l) using the apogee/perigee constraint c*p = 0 results in 
the required horizontal and radial components of velocity for apogee/ perigee 
designation maneuvers. 

= [r c 2 /(r - cos 8) ] 1 ^ 2 
(1 - l/r)/tan (8/2) 

The derivation of the equation 

0 = cos -1 [ r p cos (e^)/r T ] - e L ' 


where 


i - 


e L if e L s f 


L e L -irif e L >ir 


for computing the desired central angle 0 between two positions (r p , r p ) which 
satisfies the TPI constraints (e p . Ah) is discussed in Ref. 12. This equation is 
used in the Desired Position Routine. 



5. 


DETAILED FLOW DIAGRAMS 


Figures 6 and 7 contain the detailed flow diagrams of the main Shuttle rendez- 
vous targeting program and the general maneuver routine, respectively. The follow- 
ing six routines are called by these two programs. 

Iteration Routine 

This routine contains a Newton Raphson iterative driver based on numerically 
computed partials. The routine computes a new estimate of the dependent variable 
x and returns the old values of the error e and x. If the iteration counter c 
exceeds 15, a convergence switch is set equal to one. 

Coelliptic Maneuver Routine 

This routine computes a coelliptic velocity vector y N based on a target 
vehicle's state vector and a delta altitude. 

Phase Match 

This routine phase matches the target state vector to the primary state vec- 
tor. The controlling switch ( s phase ) equals tvo if the leading vehicle leads the 
other vehicle by more than one revolution: otherwise the switch equals one. If the 
primary vehicle leads to target vehicle, the switch is negative. 

Desired Position Routine 

This routine updates a specified state vector to the time tp and then offsets 
the updated state vector through either (AO , Ah) or (e^, Ah), depending on the 
setting of the switch s, to obtain r^. The routine contains an iterative search to 
solve the (e^ , Ah) offset problem, where e^ is defined in Figure 2 and Ah (posi- 
tive when the target orbit is above the primary) is defined as shown below. (This 
represents the TPI geometry used in Apollo and Skylab.) 


TPI point on target orbit 
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Update Routine 


This routine updates a state vector based on the update switch 

r 


s update 


s update \ 


= 1 
= 2 
= 3 

= 4 

= 5 


v. 


Updates through the time tp - t 

Updates through the time A t 

Updates through revolutions 

Updates through the angle 8 

Updates to where the orbit intersects 
the line defined by .r ^ 


Search Routine 


This routine makes the following computations depending on the setting of 

the search switch s . 

search 


(>0> 
= -1 


Finds the time of the n th apsidal crossing 
and updates the state vector to that time 

Finds the time of the next perigee crossing 
and updates the state vector to that time 


s 


search 


< 


= -2 


= -3 


Updates the state vector through the time 

t „ - t and computes the elevation angle 
r 

Finds the time associated with a specified 
elevation angle and updates the state vector 
to that time 


The detailed flow charts for these routines are shown in Figures 8 to 13. 
The iterative algorithm used to determine the time associated with the elevation 
angle is described in Ref. 8. 


Each input and output variable in the routine and subroutine call statements 
can be followed by a symbol in brackets. This symbol identifies the notation for 
the corresponding variable in the desired description and flow diagrams of the 
called routine. When identical notation is used, the bracketed symbol is omitted. 
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INPUT VARIABLES 


UNIVERSAL 

CONSTANT 



IP* ^P* — T ’ — T * *’ t ’ *1* S eng ’ m 

PROGRAM For eac h i ( l<i c total number of maneuvers): 

CONSTANTS 

s astro’ S coplan ’ s direct * s exit ’ s rnan’ 

g g g q p g 

opt * outp ’ pert ’ phase’ s search’ soln’ rdes ’ tar’ 

S S S l n e 

term’ update p’ update ^ max ’ r ’ L 

t F , 6, Ah, Ah , A8, At, Av 


c = c ^ = c i = y D = 0; i =i ;i 0 =i 

p h ; P maxg max. S 

Set f according to s gn ; a T = f g^/m 


Call Search Routine 

Input: r p> Vp, r T , v t , t, t 

' s pert j* s search.’ e L. 

Output: r pi , v pi . v T1> t 

l* e LN 

JL =i 


<^" S exit i — 

EXIT 

0 



’outp . 
\^0 


> Q i 1 = unit (v pi Xr p .); i 2 = unit (v Ti X r T _ ) 

^ y P = I Pi’ -2’ y P = — Pi ’ -2’ y T = — Ti ’ ~1 


Call Phase Match Routine 


phase > 


Input: -Pi ’ -Pi' -Ti' £Ti ,t i' s pert I ,a phase 1 


Output: r Ti , v Ti 
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At i t S + 2 

Call GMR (see below) 

°1 = Av t- At i = l s 

Call GMR (see below) 
e = (c^-Avrp)/4 



■ 




Call Iteration Routine 


Input: c, e, tg (x) , e Q , tg Q (x Q ), 4 

Output : c, t s (x), e 0 ,t s0 (x 0 ),s faU 



Call General Maneuver Routine (GMR) 


Pi’ 

-Pi’ -Ti’ 

Fi ' 

A 

r 

m, e L , a, 

1 S 1 

direct/ eng’ 

soln 

, i’ tar. 

update^ 

r 

P ( i 

+ 1) L-PN . 


— T ( i + 1 ) [— TN ] * — T ( i + 1 ) L — TN J ' 

— lc ' ±N’ A ^i’ t i+l[ t N]* Av i’ 
Av T’ e h’ e p’ S proj ' Ar proj 


Figure 6b. Main Program - Detailed Flow Diagram 
















Figure 6c. Main Program - Detailed Flow Diagram 
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V 



Figure 6d. Main Program - Detailed Flow Diagram 
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UNIVERSAL 

CONSTANTS 


INPUT VARIABLES 



(Figure 7c) 

Figure 7a. General Maneuver Routine - Detailed Flow Diagram 
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1 



6- e - sign 

(c p X r T ) * <r p X v p ) j 

-1 

cos 

[— P ‘ -T /r p r T] 


I 


Call Update Routine 


Input: r T [r], y T [v], t, t p , At, e . 
Output: r TF [r F ], Y TF [v p ] , t F 


s s 
pert' ‘ update^ 



Figure 7b. General Maneuver Routine - Detailed Flow Diagram 











— »J i = unit(v TN x r TN ); Av =-v pN • i; Av - AvT 


Figure 7c. General Maneuver Routine - Detailed Flow Diagram 
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Call Update Routine 


Input: r TN [i:] * [x] • i: pN [z D ] . *!*[*] » s pert' S update 5 

Output: [Xp] » ipp [-pi 



Figure 7d. General Maneuver Routine - Detailed Flow Diagram 
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INPUT VARIABLES 



Figure 8. Iteration Routine - Detailed Flow Diagram 


UNIVERSAL 

CONSTANTS INPUT VARIABLES 



Figure 9. Coelliptic Maneuver Routine - Detailed Flow Diagram 
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INPUT VARIABLES 



Figure 10. Phase Match Routine - Detailed Flow Diagram 
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PROGRAM CONSTANTS 


INPUT VARIABLES 



Figure 11. Desired Position Routine - Detailed Flow Diagram 
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UNIVERSAL 

CONSTANTS 


PROGRAM 

CONSTANTS 


INPUT VARIABLES 
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Figure 12b. Update Routine - Detailed Flow Diagram 
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UNIVERSAL PROGRAM 

CONSTANTS CONSTANTS INPUT VARIABLES 



(Figure ,13a. Search Routine - Detailed Flow Diagram 


5-16 

















4 



5-17 













Figure 13c. Search Routine - Detailed Flow Diagram 
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6 . 


SUPPLEMENTARY INFORMATION 


The Shuttle rendezvous targeting program proposed herein uses the basic 
targeting philosophy employed in Apollo and Skylab. (See Ref. 3 and 4. ) The 
Shuttle program represents, to some degree, a general solution to the rendezvous . 
targeting problem. The program has the capability of solving the Apollo and Skylab 
rendezvous configurations as well as many other configurations that can be deter- 
mined by specifying sets of secondary constraints. The constraints contained in 
this document may require modification as experience is gained in generating sets 
of secondary constraints based on a variety of Shuttle missions. Reference 9 
contains a more detailed discussion of the nature and number of the secondary con- 
straints employed in the maneuver sequences. Reference 11 is a users' guide to 
the Shuttle targeting program. 

The targeting program logic has been verified by a computer simulation that 
disclosed no iterative convergence problems for the range of trajectories considered. 
References 10 and 13 contain simulations of the Skylab five-maneuver rendezvous 
configuration and the six maneuver (through TPI)"Mission 2" shuttle rendezvous, 
respectively. This program is currently being integrated with the Orbiter naviga- 
tion programs to verify their combined effectiveness on a number of Orbiter mis- 
sions. 

Some of the maneuvers described in this document were based on conic orbit 
assumptions although the conic/ precision switch was set to indicate that oblate- 
ness should be considered. One example is the computation of the "circular orbit" 
maneuver based on the conic circular velocity vector. This calculation can be 
modified to provide a maneuver that would result in a minimum altitude change 
orbit in an oblate gravity field. (See Ref. 2. ) 

As TPI is generally desired at the midpoint of darkness, a stored value for 
the sun's position vector could be used to generate a display giving the angular 
distance between TPI and the midpoint of darkness. This vector could also be used 
to directly define the location of the TPI point (using the Update Routine option of 
updating to be colinear with a specified vector). 
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A decision has not been reached on a final listing of the desired astronaut dis- 
plays. As no program alarms are set for low perigees (following Skylab practice), 
the perigee altitudes for each maneuver segment contained in the maneuver sequence 
can be displayed. (If the orbits and maneuver magnitudes are approximately nominal, 
the astronaut can probably assume that the perigee altitudes will be acceptable.) 

Some standard orbit shaping maneuvers (e. g. circularization, altitude change) 
are included as rendezvous targeting maneuvers in this document. Hence, the 
rendezvous targeting program as herein contained has a limited capability to perform 
orbit shaping maneuvers assuming a dummy target vehicle state vector is used. 

The program can be expanded to allow the computation of the orbit shaping maneuvers 
for the shuttle mission by eliminating the need for the dummy state vector and adding 
additional options for orbit shaping maneuvers. However, the additional program 
complexity and the large number of inputs that have to be set nominally equal to zero 
do not seem to justify a unified approach to the calculation of the rendezvous targeting 
and orbit shaping maneuvers. Once all the orbit shaping maneuvers are defined, a 
separate Space Shuttle GN&C Equation Document will be issued for orbit shaping 
targeting. This document might utilize some of the routines herein contained. 


6-2 



REFERENCES 


1. Brand, T. , "Precision Required Velocity Deter- 
mination," Space Shuttle GN&C Equation Docu- 
ment No. 13, Rev. 1, Draper Laboratory, 

September 1971. 

2. Geyling, F. , and Westerman, H. , Introduction 
to Orbital Mechanics (Addison Wesley Co, , 

Reading, Mass., 1971), p. 212. 

3. "Guidance System Operations Plan for Manned 
CM Earth Orbital and Lunar Missions Using Pro- 
gram Colossus, " R-577, Section 5 Guidance Equa- 
tions (Rev. 14), Draper Laboratory, March 1971 

4. "Guidance System Operations Plan for Manned 
CSM Earth Orbital Missions Using Program 
Skylark I, " R-693, Section 5 Guidance Equations, 

Draper Laboratory, October 1971. 

5. Brand, T. , Brown, D. , Higgins, J. , "Unified Powered 
Flight Guidance", Space Shuttle GN&C Equation Document 
No. 24, Draper Laboratory, June 1973. 

6. Shepperd, S. , Robertson, W. , "State Extrapola- 
tion", Space Shuttle GN&C Equation Document 
No. 25, Draper Laboratory, May 1973. 

7. Robertson, W. , "Precision State and Filter Weighting 
Matrix Extrapolation, " Space Shuttle GN&C Equation 
Document No. 3, Rev., Draper Laboratory, Oct. 1972. 

8. Tempelman, W. , "The Iterative Algorithm Used in the 
Elevation Angle Search in the Rendezvous Targeting 
Program, " STS Memo No. 8-72, Draper Laboratory 
February 11, 1972. 


R-l 



9 . 


Tempelman, W. , "The Secondary Constraints Employed 
in the Orbiter Rendezvous Targeting Program, " STS 
Memo No. 52-71, Draper Laboratory, October 15, 1971. 

10. Tempelman, W. , "A Skylab Rendezvous Trajectory 
Computed with the Orbiter Targeting Program, " STS 
Memo No. 55 -71, Draper Laboratory, October 12, 1971. 

11. Tempelman, W. , "Users' Guide to the Shuttle Rendezvous 
Targeting Program., 11 STS Memo No. 9-72, Draper 
Laboratory, 14 March 1972. 

12. Tempelman, W. , "The Calculation of the TPI Central 
Angle, " Skylark Memo No. 4, Draper Laboratory, 

January 4, 1971. 

13. Tempelman, W. , "The Computation of the Mission 2 
Rendezvous Profile Using the Shuttle Rendezvous Targeting 
Program", STS Memo No. 35-73, Draper Laboratory, 

May 2, 1973. 


R-2 



